home *** CD-ROM | disk | FTP | other *** search
/ IRIX 6.2 Development Libraries / SGI IRIX 6.2 Development Libraries.iso / dist / complib.idb / usr / share / catman / p_man / cat3 / complib / slacon.z / slacon
Text File  |  1996-03-14  |  2KB  |  67 lines
  1.  
  2.  
  3.  
  4. SSSSLLLLAAAACCCCOOOONNNN((((3333FFFF))))                                                          SSSSLLLLAAAACCCCOOOONNNN((((3333FFFF))))
  5.  
  6.  
  7.  
  8. NNNNAAAAMMMMEEEE
  9.      SLACON - estimate the 1-norm of a square, real matrix A
  10.  
  11. SSSSYYYYNNNNOOOOPPPPSSSSIIIISSSS
  12.      SUBROUTINE SLACON( N, V, X, ISGN, EST, KASE )
  13.  
  14.          INTEGER        KASE, N
  15.  
  16.          REAL           EST
  17.  
  18.          INTEGER        ISGN( * )
  19.  
  20.          REAL           V( * ), X( * )
  21.  
  22. PPPPUUUURRRRPPPPOOOOSSSSEEEE
  23.      SLACON estimates the 1-norm of a square, real matrix A.  Reverse
  24.      communication is used for evaluating matrix-vector products.
  25.  
  26.  
  27. AAAARRRRGGGGUUUUMMMMEEEENNNNTTTTSSSS
  28.      N      (input) INTEGER
  29.             The order of the matrix.  N >= 1.
  30.  
  31.      V      (workspace) REAL array, dimension (N)
  32.             On the final return, V = A*W,  where  EST = norm(V)/norm(W) (W is
  33.             not returned).
  34.  
  35.      X      (input/output) REAL array, dimension (N)
  36.             On an intermediate return, X should be overwritten by A * X,   if
  37.             KASE=1, A' * X,  if KASE=2, and SLACON must be re-called with all
  38.             the other parameters unchanged.
  39.  
  40.      ISGN   (workspace) INTEGER array, dimension (N)
  41.  
  42.      EST    (output) REAL
  43.             An estimate (a lower bound) for norm(A).
  44.  
  45.      KASE   (input/output) INTEGER
  46.             On the initial call to SLACON, KASE should be 0.  On an
  47.             intermediate return, KASE will be 1 or 2, indicating whether X
  48.             should be overwritten by A * X  or A' * X.  On the final return
  49.             from SLACON, KASE will again be 0.
  50.  
  51. FFFFUUUURRRRTTTTHHHHEEEERRRR DDDDEEEETTTTAAAAIIIILLLLSSSS
  52.      Contributed by Nick Higham, University of Manchester.
  53.      Originally named SONEST, dated March 16, 1988.
  54.  
  55.      Reference: N.J. Higham, "FORTRAN codes for estimating the one-norm of a
  56.      real or complex matrix, with applications to condition estimation", ACM
  57.      Trans. Math. Soft., vol. 14, no. 4, pp. 381-396, December 1988.
  58.  
  59.  
  60.  
  61.  
  62.  
  63.                                                                         PPPPaaaaggggeeee 1111
  64.  
  65.  
  66.  
  67.